Double vector meson production from the BFKL equation 
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The double vector meson production in two photon collisions is addressed assuming that the 
color singlet t-channel exchange carries large momentum transfer. We consider the non-forward 
solution of the BFKL equation at high energy and large momentum transfer and estimate the total 
and differential cross section for the process 77 — > V1V2, where Vi and V2 can be any two vector 
mesons {Vi = p,uj,(j),J/'i',T). A comparison between our predictions and previous theoretical 
results obtained at Born level or assuming the Pomeron-exchange factorization relations is presented. 

• Our results demonstrate that the BFKL dynamics implies an enhancement of the cross sections. 

' Predictions for the future linear colliders (TESLA, CLIC and ILC) are given. 

^ . PACS numbers: 12.38.Aw, 13.85.Lg, 13.85.Ni 
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; I. INTRODUCTION 

(N ; 

^ Understanding the behavior of high energy hadron reactions from a fundamental perspective within Quantum 

Chromodynamics (QCD) is an important goal of particle physics. Attempts to test experimentally this sector of QCD 
have started some years ago with the first experimental results from lepton - hadron (DESY-HERA) and hadron - 
hadron (FERMILAB-TEVATRON) collisions. On the other hand, the papers which form the core of our knowledge 
of the Regge limit (high energy limit) of QCD were established in the late 1970s by Lipatov and collaborators 
The physical effect that they describe is often referred to as the QCD Pomeron, with the evolution described by the 
BFKL equation. The simplest process where this equation applies is the high energy scattering between two heavy 
quark-antiquark states, i.e. the onium-onium scattering. For a sufficiently heavy onium state, high energy scattering 
O ■ is a perturbative process since the onium radius gives the essential scale at which the running coupling is evaluated. 
I ' This process was proposed as a gedanken experiment to investigate the high energy regime of QCD in Refs. |^ Q 
(see also Refs. 0). At leading order (LO), the cross section grows rapidly with the energy. 
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because the number of dipoles in the light cone wave function grows rapidly with the energy. Some examples of 
■ processes which are directly related with the onium - onium scattering are forward jets in deeply inelastic events at 
low values of the Bjorken variable x in lepton-hadron scattering, jet production at large rapidity separations in hadron- 
hadron collisions and off-shell photon scattering at high energy in e~ colliders, where the photons are produced 
from the leptons beams by bremsstrahlung (For a review of BFKL searches, see e.g. Ref. 0). This last process 
presents some theoretical advantages as a probe of QCD Pomeron dynamics compared to the other ones because it 
does not involve a non-perturbative targ et 0) IE 0- In particular, the photon colliders offer unique possibility to 
probe QCD in its high energy limit The simplicity of the initial state and the possibility of study of many 

different combinations of final states making this process very useful for studying the QCD dynamics in the limit of 
high center-of-mass energy ^/s and fixed momentum transfer t. For instance, if we consider quark production in 77 
collisions we have that in the general case there can be three large momentum scales in this process - the photon 
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virtuality Q^, the momentum transfer t and the quark mass M. It imphes that for real photons interactions and 
t = only heavy quark production can be treated perturbatively. On the other hand, perturbation theory can be 
used even for light quark production if \t\ 3> ^qcd- ^^^^ years many authors have studied in detail the heavy 

quark production in 77 collisions considering different theoretical approaches (See also Ref. 

Other possibility for the study of the QCD Pomeron is the vector meson pairs production in 77 collisions • At 
very high energies s 3> — diffractive processes such as 77 neutral vector (or pseudo-scalar) meson pairs with 
virtual or real photons can test the QCD Pomeron (Odderon) in a detailed way utilizing the simplest possible initial 
state. As in the case of the large angle exclusive 77 processes, the scattering amplitude is computed by convoluting 
the hard scattering pQCD amplitude for 77 —^ qqqq with the vector meson wave functions. For heavy vector mesons, 
this cross section can be calculated using the perturbative QCD methods. First calculations considering the Born 
two-gluon approximation have been done in Refs. 1 141 . Some years ago, the double J/^* production in 77 collisions has 
been proposed as a probe of the QCD Pomeron [l^, where the cross sections are calculated in the BFKL framework 
and also considering a next leading order (NLO) calculation by a kinematic veto on the gluon momenta in the 
ladder. Non-perturbative contributions for this process are studied in |16j by employing modified gluon propagators. 
Theoretical estimates of the cross sections presented in llSl have demonstrated that measurement of this reaction 
at a Photon Collider should be feasible (See also Ref. On the other hand, in order to calculate light vector 

meson production it is necessary to consider virtual photon collisions or a large momentum transfer. In this paper 
we study the double vector meson production in 77 collisions, where the the color singlet t-channel exchange carries 
large momentum transfer. Such processes are characterized by a large rapidity gap between the vector mesons, 
with it identified experimentally by their decay into charged particles, e.g. p — > tt+tt". The presence of a large 
momentum transfer allow us to calculate the light and heavy vector mesons production cross sections for real photon 
interactions. We consider the non- forward solution of the BFKL equation at high energy and large momentum transfer 
and estimate the total and differential cross section for the process 77 V1V2, where Vi and V2 can be any two vector 
mesons (Vi = p, w, 0, J/^', T). Our study is motivated by the fact that the experimental data for the vector meson 
photo-production at high t in electron-proton collisions at HERA can be described quite well using the impact factor 
representation and the non- forward BFKL solution [T9I l20j . In order to simplify our calculations we will consider the 
non-relativistic approximation in the calculations of the impact factors. This approximation should be not applicable 
for light mesons. However, as verified in Ref. jl9j . the predictions obtained using the non-forward BFKL solution at 
asymptotic energies (leading conformal spin solution) and this approximation agree quite well with heavy and light 
meson HERA data. This phenomenological success occurs due to the selection of the dominant amplitude of the 
general expression as demonstrated in Ref. 20]. We will assume that this dominance will be also present in the 
processes considered in this paper, and use this approximation in that follows. W"e postponed for a future publication 
the study of double vector production using the full amplitude, similarly to Ref. '20'| , which will allow us consider the 
helicity flip of the quarks. 

This paper is organized as follows. In next section we present a brief review of the formalism necessary in order to 
calculate the differential and total cross section of the 77 V1V2 process. In Section lHll we present our predictions and 
compare our results with previous estimatives. Finally, in Section llVI we present a summary of our main conclusions. 



II. FORMALISM 



In order to calculate the 77 ^ V1V2 cross section we will consider the impact factor representation proposed by 
Cheng and Wu many years ago, which allow us factorize the process in three disjoint parts: the impact factors 
associated to the transitions X^^' (i = 1,2) and the exchange mechanism between them represented by JCbfkl in 
Fig. n In this representation, the amplitude for the high energy process AB —> CD can be expressed on the form 

A^^-^^iq) = / d'kd'k' I^-^ik,q) I^^^(fe',q), (1) 

J k {q~ k)^ 

where X^^'~^ and I^^^ are the impact factors for the upper and lower parts of the diagram, respectively. That is, 
they are the impact factors for the processes A C and B ^ D with two gluons carrying transverse momenta k 
and q — k attached, the gluons being in an overall color singlet state. At lowest order the process is described by two 
gluon exchange, which implies JCbfkl 5^'^\k ~ k') and an energy independent cross section. At higher order, the 
dominant contribution is given by the QCD pomeron singularity which is generated by the ladder diagrams with the 
(reggeized) gluon exchange along the ladder. The QCD pomeron is described by the BFKL equation [ij, which implies 
that the exchange of the gluon ladder with interacting gluons generates increasing cross sections with the energy. As 
our goal is the analysis of double vector meson production at large - t, we will use in our calculations the non-forward 
solution of the BFKL equation in the leading logarithmic approximation, obtained by Lipatov in Ref. ^3|- 



The differential cross-section, characterized by the invariant coUision energy squared s, is expressed in terms of the 
amphtude A{s,t) as foUows 

The amphtude is dominated by its imaginary part, which we shaU parametrize, as in |l9l l2^ . by a dimensionless 
quantity 

ImA{s,t)^^n^,T) (3) 



where z and r are defined by 



where My is the mass of the vector meson, is the photon virtuality and is a characteristic scale related to My 
and \t\. In this paper we only consider Qj = 0. In LLA, A is arbitrary (but must depend on the scale in the problem, 
see discussion below) and as is a constant. For completeness, we give the cross-section expressed in terms of J-{z,t), 
where the real part of the amplitude is neglected, 

da (jyy ^ V1V2) l^"" IT-/ m2 



dt 81i4 
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This representation is rather convenient for the calculations performed below. 

The BFKL amplitude, in the leading logarithm approximation (LLA) and lowest conformal spin (n = 0), is given 
by M 

^BFKL(z,r) = d.j-^^^e^^^^^r^^{Q^)irHQ±r, (7) 

where Q± is the momentum transfered, t — —Q\, (the subscript denotes two dimensional transverse vectors) and 

x(^^)=47^c(^^(l)-V^Q+^^.^^ (8) 

is proportional to the BKFL kernel eigenvalues p^ . with "0(2;) being the digamma function. 

The quantities /J^' are given in terms of the impact factors X^y. and the BFKL eigenfunctions as follows [2^ . 

/r'(Q±) = J ^i^v.{k^,Q±) I d^pid^P2 (9) 
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In the case of coupling to a colorless state only the first term in the square bracket survives since XA{k±, Q± — k±) — 
lA{k± = 0,Q±) = in this case. The impact factor T~^Vi describes in the high energy limit the couplings of the 
external particle pair to the color singlet gluonic ladder. It are obtained in the perturbative QCD framework and we 
approximate them by the leading terms in the perturbative expansion psj : 

Citts / 1 



In this formula, factorization of the scattering process and the meson formation is assumed, and the non-relativistic 
approximation of the meson wave function is used. In this approximation the quarks in the meson have coUinear four- 
momenta and Alv — 'i.Mq where Mq is the mass of the constituent quark. To leading order accuracy, the constant Ci 
may be related to the vector meson leptonic decay width 

= (11) 



Moreover, we have that 



= ql + Qi/4, (12) 
ql = (Q2+m2.)/4. (13) 

Ahhough this approximation should be not apphcable for Hght mesons we will use it here following p^ . which has 
observed that this approximation implies a reasonable description of the heavy and light meson HERA data. As 
emphasized in Ref . "20(| , this phenomenological success occurs due to the selection of the dominant amplitude of the 
general expression. We will assume that this dominance will be also present in the processes considered in this paper. 
We postponed for a future publication the study of double vector production using the full amplitude (also considering 
higher conformal spins contributions). 

Using Eq. ^ into ®, one obtains [2II2II 

X 



Ql r(l/2 + \ A ) 7i/2_,^ 2TTI V4M2 
r2(l/2 + u)r(l/2 - u/2 - iv/2)T{l/2 - u/2 + iv/2) 
r(l/2 + u/2 - w/2)T{l/2 + u/2 + iv/2) 



The differential cross section can be directly calculated substituting the above expression in Eq. lO and evaluating 
numerically the integrals. The total cross section will be given by 

where |i|mm is the minimum momentum transfer, which we will assume different values for the distinct processes. 



III. RESULTS AND DISCUSSIONS 



Lets start the analysis of our results discussing the choice for the parameters as and A. We have that as appears 
in two different places in our calculations: in Eq. (|15|l and in the definition of the variable z [Eq. Q)]. The first factor 
comes from the couplings of the two gluons to the impact factor and the second one is generated by the gluon coupling 
inside the gluon ladder (For details see Ref. l2Qj) . In this paper we treat these strong couplings as being identical 
and assume a fixed Og, which is appropriate to the leading logarithmic accuracy. Furthermore, as the cross section is 
proportional to a^, our results are strongly dependent on the choice for the as value. We assume as = 0.2, following 
Refs. flfll. f2(Tj|. where this value was determined from a fit to the HERA data (For a detailed discussion see Ref. 19]). 
Similarly, at leading logarithmic approximation, A is arbitrary but must depend on the scale in the problem. In our 
case we have that in general A will be a function of My and/or t. Following Ref. IS] we will assume that A can 
be expressed by A^ = PMy + 7|t|, with /3 = 1 and 7 = 0. In order to consider the possibility of the production of 
two different mesons we will generalize this expression for A^ = (3iMy_^ + /?2-My^ + 7|i|. In Fig. |21we analyze the 
dependence of the differential cross section for pJ /^ production in our choice for the [3i parameters. Although we 
only present the results for 77 collisions a,t W = 500 GeV, the dependence is similar for other center of mass energies. 
As A is present in the definition of the variable z [Eq. Q], which is associated with the energy dependence of the 
cross section, we have that a different prescription for A basically modify the normalization of the differential cross 
section. In that follows we will assume /3i = /32 = 1/2. 

In Figs. 121 and ^ we present our results for the momentum transfer dependence of the differential cross sections for 
the production of different combinations of pairs of vector mesons and distinct values of energy. For comparison, the 
Born level prediction which corresponds a two gluon exchange, given by the following formula, 

^-""•"-^"VlFMlf^ (16) 

is also presented. We have that the differential cross sections decreases when the masses of the pair of mesons increases 
as well as it becomes flatter. Another remarkable feature is the change of the behavior of the cross section when we 
take into account the BFKL machinery in comparison with the Born level cross section. The BFKL dynamics implies 
a growth of the differential cross section and a different slope in comparison to Born level, which is not modified 
when the energy increases. As already pointed out in Ref. [13, the inclusion of BFKL dynamics result in a steeper 
i-dependence than the Born calculation. 



The total cross section for different combinations of mesons are obtained directly from Eq. H15|l and the results are 
displayed in the Figs. El El and and in the Tables HI and ITU For light vector meson production, a lower cut-off in the 
t-integration is necessary in order to minimize non-perturbative (Soft Pomeron) contributions. We choose in this case 
\t\min — 1 GeV^, which is reasonable considering that the HERA data for the photoproduction of light vector mesons 
in this kinematic region are quite well described using a similar approach |l9l l20j|. On the other hand, we assume 
l^lmm = when at least a heavy meson is produced, since in this case we have a hard scale present which justifies 
the perturbative calculations in the \ow-t {t ~ 0) region. The predictions obtained using \t\„iin — 1 GeV^ are also 
presented in Figs. El and El We have that the total cross sections are smaller for heavier vector mesons. Moreover, 
the introduction of a cut-off in the f-integration implies a reduction of the normalization of the total cross section, 
as expected from the analyses of the differential cross section. Another important feature is that difference between 
the calculations assuming |t|mm = and |i|mm = 1 GeV^ is also smaller for heavier vector mesons. In Fig. [3 we 
present our predictions for the production of two light vector mesons assuming |i|„im — 1 GeV^. We have that all 
cross sections present a strong growth with the energy [atot oc W^), with a similar slope (A « 1.4). 

In the Tables HI and Hll we present a comparison between our predictions for the total cross section for some values 
of center of mass energy and the energy independent prediction which comes from the Born amplitude. We have that 
the inclusion of the BFKL dynamics implies that the cross sections strongly increases with the energy, resulting in an 
enhancement of the cross section in comparison with the Born level prediction. This enhancement is larger when a 
heavy vector meson is produced and \t\min = GeV^. If we assume \t\min = 1 GeV^ the BFKL dynamics predicts a 
similar enhancement of the light and heavy vector meson cross sections. 

Lets now compare our results with those obtained in other approaches 0, 0, 0, 0, |23, 113, |^ (For related 
discussions see Refs. 113, US, 113 the FELIX collaboration has proposed the construction of a full 

acceptance detector for the LHC, with a primary proposal providing comprehensive observation of a very broad range 
of strong-interaction processes. In particular, studies of two photon physics in AA collisions has been discussed in its 
proposal ^2l\. There, the differential cross-section da{-fj ViV2)/dt, with Vi = p, J was parameterized in the 
form 

— (77^FiF2) = v4y,y, ( — j expUBv,v,+4ta'y^yJn—\. (17) 

This parameterization was obtained assuming the validity of the pomeron-exchange factorization at high energies 
(See also Ref. jll]) which relates the cross sections for different processes. The parameters AviV2j By^Vi jCviVi and 
a'y^ V2 were fixed considering the data for vector meson photoproduction at HERA and for pp elastic scattering. For 
the pJ/'^ case, for instance, Apjj^ ^ 1.1 nb/GeV^^, Epj/^, ~ 2.5 GeV^^, Cpjj^ = 0.8 and a^j^^ — 0. We have 

assume this parameterization and calculated the total cross section considering \t\m.in — and \t\min = 1 GeV^ (See 
Fig. EJ. In comparison with the BFKL predictions we have that the FELIX parameterization implies a different 
normalization and energy dependence, with the BFKL one being steeper. Basically, the FELIX parameterization 
overestimate the cross sections at low energies and underestimate the cross section at very high energies. The only 
exception occurs when we consider the double J j'^ production. In this case the BFKL prediction is larger than the 
FELIX one for all energies if we assume Ajj^jj^ = 3.1 x 10"'' nb/GeV~^, Bj/q,j/q, ~ 1.5 GeV~^, Cjj^jj^ = 1.38 
and a'j/^j/^ = 0. On the other hand, if we assume Aj/^j/^ = 0.16/ \n{W/Mj/^) nb/GeV~^ and Cj/^j,j/^;, — 2.0 as 
also proposed in Ref. [131 (BFKL-fit) we have that the FELIX parameterization implies a larger cross section for all 
energies. The main conclusion of the comparison between the BFKL predictions and those obtained assuming the 
pomeron-exchange factorization is that the BFKL dynamics implies a violation of this factorization, which could be 
tested in future colliders. 

In Ref. [1^ the double vector meson production in 77 collisions has been analyzed using the stochastic vacuum 
model |35| | . This model is based on the assumption that the infrared behavior of QCD can be approximated by 
a Gaussian stochastic process in the gluon field strength tensor. A shortcoming of this model is that the energy 
dependence is not predicted, which implies the introduction of a p henomenological ansatz for this dependence. In 
Ref. [isl the authors have considered the two-Pomeron model j3(il |. resulting that the scattering amplitude is given 
by the general expression 

where asoft{t) = 1-08 -I- 0.25 1 and ahard{t) = 1.28. These values are used in order to describe the pp and ep data. It is 
important to emphasize that as this model includes soft and hard contributions the light and heavy vector meson total 
cross sections can be calculated assuming \t\,nin — GeV^. Consequently, its results for double light vector meson 
production cannot be directly compared with our predictions. However, analyzing only the energy dependence of 
both predictions, we have that the BFKL dynamics implies a steeper dependence as ex pec ted, since the soft Pomeron 
contribution dominates the cross sections for these processes in the model proposed in [l8|. On the other hand, when 



a heavy vector meson is present a comparison can be made. In this case the predictions agree for low values of energy 
{W w 10 GeV) but are distinct for larger energies due to the larger intercept resulting from the BFKL equation than 
the phenomenological one ahardit)- 

The pJ/'i production in 77 processes has also been estimated in Refs. |30ll3ll |. There, the differential cross section 
was estimated in a similar way as the elastic J/^* photoproduction off the proton [2^ and reads as [sof. 

^.h-,-^,Jm ,c '^-'I'-l'^y ' 1.G>(..MJ,,/4)1^ exp(B„„.<) , (19) 

where C denotes factors of corrections discussed in detail in Refs. [sill l8ll | and x — Mj^^/W^^. In the small-i 
approximation, the slope is estimated to be Bp j/^ — 5.5 ± 1.0 GeV^^. The light meson-photon coupling is denoted 
by = 0.454. The process above was proposed as a probe of the gluon distribution on the meson xC^ and, by 
consequence, a constrain for the photon structure. The enhancement in the sensitivity by taking the square of 
those distributions in the total cross section could discriminate them in measurements at the future photon colliders. 
Considering different models for the gluon distribution of the meson the total cross section was estimated as being 
almost 10 pb at « 10 GeV and in the range between 100 and 1000 phaXW^i 100 GeV. Our results agree with these 
predictions, with a behavior similar to those obtained using the GRS(LO) parameterization for the gluon distribution 
on the p meson (For details see Refs. [s^, 0|)- Moreover, in Ref. 31] the double p production was estimated 
assuming the pomeron-exchange factorization to relate distinct cross sections, similarly to FELIX collaboration. Our 
results agree with its results. However, as its predictions are strongly dependent on the assumptions present in the 
calculations of the double J/'^ and pJ/"^ production (See Table 1 in Ref. ^^), a direct comparison is not very 
illuminating. 

The first computations of the double vector meson production at arbitrary momentum transfer were done in Refs. 
[T^ considering the Born graph for the Pomeron (two- gluon exchange). The total cross section for double light 
meson production were computed assuming \t\min = 3 GeV^ and as = 0.32. Taking into account these differences 
we have that our calculations of the Born level reasonably agree with those presented in Refs. 0|. The double J/vP 
production has been estimated in Refs. [l^lisllr^ at Born level and assuming the BFKL dynamics. Similarly to light 
vector meson production, our results agree with the predictions from [l4| if the differences are taken into account. 
In Ref. the total cross section was estimated assuming a small-t approximation and an ansatz for the -Bj/^j/^ 
slope parameter. Moreover, it was calculated considering the solution of the forward BFKL solution and some of its 
next-to-leading corrections. The normalization of the cross section is strongly dependence of the ag value (See Table 
1 in Ref. IJj) and its energy dependence is steeper when the leading order (LO) solution of the BFKL equation is 
used. Basically, although the normalization of the distinct predictions differ, the energy dependence is very similar. 
The difference in the normalization can be explained in terms of the uncertainty introduced in the calculations by the 
small-t approximation used in Ref. [17ij . On the other hand, if we compare our results with those obtained using the 
solution of the forward BFKL equation with a different intercept, we have that the normalization for low energies is 
similar, but the energy dependence of those predictions are softer. This result is expected, since the next-to-leading 
corrections for the BFKL equation implies a lower intercept than the leading order solution used in our calculations. 

Finally, lets estimate the expected number of events of some of the processes calculated in this paper for the future 
linear colliders. Currently, the international High Energy Physics (HEP) community is considering the construction 
of a new e+e~ collider |33,|3l,|33- For 77 collisions with center of mass energies equal to W = 500 GeV, luminosities 
of order C = 1.1, 0.2 and 4.0 x 10^"* cm~^s~^ are expected at TESLA, CLIC and ILC, respectively. 

It is important to emphasize that the next generation of linear collider can reach values of center of mass energies 
up to 3000 GeV. As the BFKL dynamics predicts a strong growth of the cross sections with the energy, the number 
of events presented in that follows may be considered a lower bound. In Table Ulll we present our predictions for the 
number of events per second for the pp, pJ/"^, J/'i'J/'^ and TT production. For double p production our estimate 
is conservative, since the soft Pomeron contribution is not included in our calculations. We predict a large number 
of events related to double meson production in 77 collisions, allowing future experimental analyses, even if the 
acceptance for vector meson detection were low. Consequently, we believe that this process could be used to constrain 
the QCD dynamics at high energies. 



IV. SUMMARY 



The photon colliders offer unique possibility to probe QCD in a new and hitherto unexplorated regime. Moreover, 
good knowledge and understanding of two photon processes will be essential for controlling background contributions 
to other processes. In this paper we have studied the double vector meson production using the non-forward solution 



of the BFKL equation at high energy and large momentum transfer. The total and differential cross section for 
the process 77 — > V1V2, where Vi and V2 can be any two vector mesons (Vi = p, w, 0, J/^P, T), were estimated 
and compared with those obtained using different approaches. In particular, we have verified that the pomeron- 
exchanged factorization leads to different predictions than those obtained in this paper. Our main conclusion is that 
the forthcoming photon colliders could experimentally to check our predictions. We believe that this process could 
be used to constrain the QCD dynamics at high energies. However, several points deserves more detailed studies: the 
contribution of all conformal spins, the helicity flip of the quarks, the next-to-leading order corrections to the BFKL 
dynamics and the corrections associated to the saturation effects. We plan to estimate these contributions in future 
publications. 
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Born 


W = 100 GcV 


W = 500 GcV 


W = 1000 GeV 


pp 


3.18 


51.0 


730.0 


2019.0 


fJLO 


0.27 


4.39 


59.0 


173.0 


p(j) 


0.48 


9.52 


135.0 


367.0 


LOU) 


0.023 


0.37 


4.9 


14.0 


LU<f) 


0.04 


0.81 


11.6 


31.0 


H 


0.076 


1.87 


25.2 


71.0 



TABLE I: The double vector meson production cross sections in processes at different energies and |i|m.m — 1 GeV^ Cross 
sections are given in pb. 





Born 


W = 100 GeV 


W = 500 GeV 


W = 1000 GeV 




1.43 (0.29) 


750.0 (6.8) 


2550.0 (89.0) 


4833.0 (234.0) 




0.12 (0.02) 


24.0 (0.59) 


410.0 (5.6) 


1199.0 (20.0) 




0.21 (0.05) 


34.0 (1.6) 


460.0 (17.0) 


1647.0 (55.0) 


J/*J/* 


0.19 (0.10) 


13.0 (3.6) 


190.0 (42.0) 


586.0 (125.0) 


pT 


0.005 (0.001) 


0.47 (0.03) 


6.7 (0.35) 


24.0 (1.1) 


wT 


0.0005 (0.0001) 


0.04 (0.003) 


0.75 (0.03) 


2.09 (0.1) 


4,T 


0.001 (0.0003) 


0.06 (0.009) 


1.1 (0.09) 


3.25 (0.29) 




0.0025 (0.0019) 


0.07 (0.03) 


0.88 (0.45) 


2.91 (1.27) 




D.OOOl:! (().()()()12) 


I). 002 (O.OOi) 


0.02 (O.OLS) 


O.OG (0.05) 



TABLE II: The double vector meson production cross sections in 77 processes at different energies assuming \t\min = 0(1) 
GeV^. Cross sections are given in pb. 





TESLA 


CLIC 


ILC 


PP 


8.0 


2.0 


28.0 


p.7/* 


28.0 


5.0 


100.0 


J/*,//* 


2.0 


0.1 


8.0 




0.0002 


0.00004 


0.0007 



TABLE III: Number of events per second for double vector meson production at TESLA, CLIC and ILC luminosities (W = 
500 GeV). For the double p production we assume \t\min = 1 GeV^. 




FIG. 2: Dependence of the differential cross section for pJ production in the choice for the /3 parameters. 
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FIG. 3: Momentum transfer dependence of the differential cross sections for the production of different combinations of pairs of 
vector mesons and distinct values of center of mass energy. Sohd line: W = 100 GeV; dashed line: W = 500 GeV; dot-dashed 
line: W = 1000 GeV; dot-dot-dashed: Born level. 




FIG. 4: Momentum transfer dependence of the differential cross sections for the production of different combinations of pairs of 
vector mesons and distinct values of center of mass energy. Solid line: W = 100 GeV; dashed line: W = 500 GeV; dot-dashed 
hue: W = 1000 GeV; dot-dot-dashed: Born level. 
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FIG. 5: Total cross section for the production of different combinations of pairs of vector mesons. Solid line: BFKL (|t|min = 
GeV^); dashed line: BFKL {\t\min = 1 GeV^); dot-dashed line: FELIX {\t\min = GeV^); dot-dot-dashed line: FELIX 
(|t|mm ~ 1 GeV^). The dot-dashed-dashed line corresponds to the BFKL-fit from FELIX (see text). 




